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Supersonic/hypersonic laminar heating correlations that were developed for damage- 
assessment analysis of atmospheric re-entry vehicles have been modified and extended to 
cover fully-turbulent conditions over rectangular cavity geometries that are aligned with the 
local velocity. Turbulent boundary layer properties were computationally determined and 
used to develop the cavity geometry parametrics and to correlate experimental closed cavity 
heating data to yield new relationships for the floor-averaged and centerline endwall peak- 
heating augmentation. With the form of the closed-cavity correlations established, historical 
data were used to develop new correlations for turbulent open-cavity heating. 


Nomenclature 

a, b correlation curve fit coefficients, see Eq. 2 
BF Bump Factor 

h heat transfer coefficient, h=q/(H aw -H w ), (lbm/ft 2 /s) 

F, W, H cavity length, width, and depth (in) 

FCF/UCF lower/upper confidence limit 

FPF/UPF lower/upper prediction limit 

M Mach number 

M e boundary layer edge Mach number at the cavity entrance 

R 2 Finear Correlation Coefficient 

Re unit Reynolds number (1/ft) 

Re e momentum thickness Reynolds number 

X, Y correlation parameters, see Eq. 2. 

x axial distance from model leading edge (in) 

y spanwise distance from model centerline (in) 

z distance normal to x-y plane (in) 

a angle-of-attack (deg) 

a, (3, a, T Bump Factor correlation coefficients, see Eq. 1 

5 boundary layer thickness (in) 

0 boundary layer momentum thickness (in) 

Subscripts 

AVG average value in cavity 

FR stagnation point conditions from Fay-Riddell calculation for a hemisphere 

FAM laminar 

MAX maximum value on cavity centerline immediately downstream of endwall 

TURB turbulent 

oo freestream static conditions 
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I. Introduction 

The August 2003 Final Report of the Space Shuttle Columbia Accident Investigation Board 1 (CAIB) 
recommended the development of improved damage assessment models for application to future flights, primarily 
for laminar flow conditions. In response, the Orbiter Entry Aeroheating Panel at the NASA Johnson Space Center 
established the Cavity Heating Team to develop the Cavity Heating Tool (CHT). This tool defines the heating 
augmentation to the vehicle outer mold lines caused by impact damage that does not penetrate the vehicle aluminum 
substructure. A survey of the cavity heating literature conducted for CHT development revealed that much general 
information was available on cavity flow physics and heating; however, information required to address Orbiter- 
specific damage was very sparse. For example, there was no information regarding the real-cavity effects of random 
geometry, cross-flow, or pressure gradient. This sparseness led to an extensive series of supersonic/hypersonic 
idealized-cavity-flow experiments by Everhart, et al. 2 " 6 that were conducted in the conventional 20-Inch Mach 6 and 
31 -Inch Mach 10 Wind Tunnels at the NASA Langley Research Center. The idealized cavity geometries were 
rectangular boxes that were easily defined by depth (H), width (W), and length (L), and scaled with the cavity- 
entrance boundary layer thickness (8). Heating was measured using global phosphor thermography methods. The 
reference documents contain an in-depth presentation of the literature survey; and, they describe the multitude of 
idealized cavity geometries tested to establish the effects of length, width, depth, entry/exit/sidewall angles, 
orientation, surface curvature, and the fluid dynamic effects of Mach number, Reynolds number, and boundary layer 
thickness. Preliminary analysis of these data enabled the development of the assessment methods that form the 
foundation of the CHT as reported by Anderson, et al. 7 

A discussion of the cavity flow physics and a more refined analysis of the laminar-entry heating augmentation or 
bump factor (BF) was presented in Reference 8 for open-flow (short) cavities where the flow skips over the cavity 
and for closed-flow (long) cavities where the flow enters the cavity and impinges on the floor. Open cavities have a 
nominal L/H<10, while closed cavities are nominally L/H>14. The transitional-cavity region between L/H=10 and 
14 may be unsteady as the flow oscillates back and forth between the two steady flow regimes. A combined single- 
curve parametric variation for the floor-averaged BF and centerline endwall maximum BF was derived, as was a 
model of the heating-wake decay downstream of the cavity. Additionally, uncertainty limits on the correlation curve 
coefficients and prediction limits about the correlation curves were also provided. 

A limited number of closed-flow turbulent test cases were obtained as a part of the Mach 6 cavity heating 
experiments to baseline heating levels and variations (Ref. 3). Using this data, the above-cited laminar correlations 
have now been extended to include turbulent boundary layer conditions. The turbulent data were originally reduced 
and reported using laminar computational solutions for entry reference conditions because turbulent solutions were 
unavailable. Turbulent solutions have since been developed, and the turbulent data have been re-reduced with these 
new entry conditions. The first task in this paper is to present a sample of the re-reduced turbulent data. Next, the 
laminar correlation parameters will be tested to determine their applicability to characterize closed-cavity turbulent 
heating variations. Coefficient values in the parameters will then be adjusted to minimize turbulent correlation 
uncertainties. Finally, using this functional relationship as a guide, correlations will be developed for turbulent open- 
flow cavities using archival data obtained from the literature. 

II. Experimental Methods 

This section presents a brief discussion of the experimental methods for the data obtained in Test 6888 (denoted 
T6888). The test facility is first discussed, followed by a description of the models and cavity design, and an 
overview of the phosphor thermography data system used to acquire the global heating measurements. Finally, the 
data reduction and presentation parameters are described. 

A. Facility 

The NASA Langley 20-Inch Mach 6 Tunnel was used to develop the turbulent cavity-heating database. This 
conventional, perfect-gas blowdown facility uses air as a test medium. It is described in detail in Micol. 9 Flow 
properties were determined using the GASPROPS code developed by Hollis. 10 

B. Models 

A schematic showing the model coordinate system orientation is presented in Figure 1. The x axis is measure 
down the length of the model and the z axis is measure normal to the surface. Initially, the origin is at the model 
leading edge; it is transferred to the cavity leading edge for data analysis. A strut-mounted 10-inch-wide by 28-inch- 
long steel flat-plate model having a blunted leading edge with a radius of 0.125 inches was used for experiments 
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(see Figure 2). The blunt nose model provides edge Mach numbers of approximately 2.9 to 3 (see Ref. 2) at zero 
incidence as required to match Orb iter flight surface conditions. The model includes an upstream insert region where 
boundary layer trips can be installed. Boundary layer trip hardware and tripping methods developed by Berry, et. 
al. 11,12 were used to ensure turbulent flow over the model. Using Berry’s specification, the height of the trip, k, is 
determined from k/8=0.6 for the present flow conditions. Off-design test conditions were examined to ensure fully 
tripped turbulent conditions. A discussion of the boundary layer trips and the baseline heating data are presented in 
Ref. 3. The model also includes a 4-inch wide by 18-inch long aft region starting 7.5 inches from the nose for test 
articles. A representative example of the T6888 cavity inserts (see Ref. 3) is shown in Figure 3. They were 
fabricated using ceramic casting methods as described in Buck, et al. 13,14 After fabrication, the cavity models were 
coated with a nominal 27 qm thick mixture of phosphors suspended in a silica-based colloidal binder 15 and sent to 
quality assurance for measurement and application of small circular locating markers, known as fiducial marks. 
Cavity size parameters and location on the insert plate are given in Table 1. 

C. Cavity Designs and Test Conditions 

Impact-induced cavities must be geometrically and fluid-dynamically scaled for damage assessment modeling. 
For the Shuttle Orb iter, a maximum cavity depth constraint can be established based on the thickness of the thermal 
protection system tiles. Because of the launch configuration, impact damage on the forward portion of the vehicle is 
likely to be deep and short, in contrast to long, shallow grazing damage that is more likely on the aft portion of the 
vehicle. Parametric variations used to establish the cavity-heating database were guided by these general 
considerations. 

Typical surface flow conditions for the windward surface of the Shuttle Orb iter were obtained from Campbell, et 
al. 16 to establish the parameter variations at the cavity entrance, (M e , Re 0 , 8) F ii g ht- Centerline test article surface 
conditions were determined for the tunnel freestream conditions from two-dimensional centerline contour 
simulations of the baseline geometry using the LAURA code. 17,18 The range of the computationally determined 
(M e , Re e )Tunnei conditions was matched to flight. The boundary layer thickness ratio (8 Tu nnei/§Fiight) was used to scale 
flight damage scenarios to tunnel scale. Cavity conditions (M e , Re 0 , 8) Tu nnei were taken from the computational 
solutions at the cavity entrance. 

All turbulent-entry data were acquired using closed (long) baseline geometries, each with a flat bottom and 
rectangular planform and cross section. Tripped turbulent flow was developed over the model surface as discussed 
in Ref. 3. Specifics of the flight-to-tunnel scaling, cavities, and their designs are included in the database 
publications by Everhart, et al., using methods and hardware developed by Berry, et al. to ensure fully-turbulent 
flow. 

D. Phosphor Thermography Technique 

Global surface heating distributions were calculated using the two-color, relative-intensity phosphor 
thermography method. 15 This is the standard technique for obtaining aeroheating data in NASA Langley’s 
hypersonic wind tunnels, and it can be used to identify the surface heating effects of complex three-dimensional 
flow phenomena, which are difficult to examine using conventional discrete-sensor methods. With this method, 
ceramic wind tunnel models are coated with phosphor crystals that fluoresce in the red and green regions of the 
visible light spectrum when illuminated by ultraviolet (UV) light. The phosphor-coated model is exposed to the 
heated flow during a wind tunnel run, and the resulting changes in fluorescence intensity of the model are recorded 
and digitized. The surface temperature distributions are determined from the fluorescence intensities through prior 
calibrations, reduced to surface heating distributions, and analyzed using the IHEAT 15 code, which implements one- 
dimensional, semi-infmite-solid heat conduction theory. These distributions may be mapped onto a 3 -dimensional 
CAD representation of the test article using the MAP3D 15 code. 

E. Data Reduction 

Global mappings of the surface heating were converted to a cavity-based reference system by translating the 
origin from the leading edge of the model to the leading edge center of each cavity. They were then scaled using the 
cavity depth and the boundary layer thickness at the cavity leading edge to yield a scaled-geometry cavity of depth 
H/8, width W/H, and length L/H. 

Local heating augmentation or bump factor (BF) was used to assess the effect of the cavity on the undisturbed 
environment. Bump Factor is defined as the local heating (hi oca i) normalized by an averaged reference heating 
(Iiavg), yielding BF=hi 0C ai/hAVG- For analysis consistency, h A vG is computed from the measured data, in a scaled- 
space region ahead of each cavity across the cavity width over the length -3 <x/H<-l. In this format, BF=1 
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becomes the nominal, fully developed, undisturbed condition for a flat plate. This location also avoids a potential 
bias caused by an expansion-induced heating peak upstream of the cavity leading edge. Two metrics, BF A vg and 
BFmax, are used to characterize the cavity heating. BF AVG was determined by averaging all imaged data inside the 
entire cavity. Only a small number of non-floor pixels were present because the upstream endwall was not visible in 
the camera imaging and only 2-3 pixels were visible in any longitudinal linecut on the downstream cavity endwall. 
Likewise, only one of the sidewalls was imaged and, typically, only 1 -2 pixels were visible in any spanwise linecut. 
This area averaging greatly reduced noise in the data caused by the low temperatures experienced on the cavity floor 
for laminar flow conditions. BF M ax was determined as the maximum centerline value; this value was typically 
located immediately downstream of the cavity endwall. 

III. Computational Simulation of Baseline Model Surface Flow 

Boundary layer edge properties were computed using the LAURA 18,19 (Langley Aero thermo dynamic Upwind 
Relaxation Algorithm) code. LAURA is a computational fluid dynamics (CFD) tool specialized for hypersonic 
reentry physics and chemistry. In the past, the LAURA code has been applied and validated for a range of vehicle 
shapes and flow conditions including the Shuttle Orbiter laminar aerothermodynamic flight data 20 , blunt planetary 
configurations, ’ and slender lifting bodies. The code is a three-dimensional upwind finite-volume solver for the 
inviscid Euler equations, thin-layer Navier-Stokes, or full Navier-Stokes governing equations. Roe-averaging 24 with 
Harten’s entropy fix 25 and Yee’s Symmetric Total Variation Diminishing (TVD) limiter 26 is used to formulate the 
inviscid fluxes, while a second order scheme solves the viscous fluxes. Point-relaxation is used to begin the iterative 
solution while a line relaxation approach is invoked to complete the integration in time to steady state. 

The Navier-Stokes equations were solved in simulating the turbulent flow over a flat-plate at ground-based test 
conditions in conjunction with a perfect gas air thermodynamic model. For this work, it was sufficient to start the 
turbulent flow at the plate leading edge and it was modeled using the Cebeci- Smith algebraic turbulence model. 27 A 
sample of the turbulent heating prediction is compared with experimental measurements for Run 31 in Figure 4a 
(M e =2.93, Re 0 =15O). Corresponding longitudinal predictions for the M e , Re 0 , and 8, are presented in Figure 4b. 
Laminar solutions are also presented in the figures for comparison. 

IV. Test Data and Parameter Space 

The turbulent-run test matrix is presented in Table 1, along with the scaled geometry and the BF AVG and BF M ax 
values for each set of conditions. Typically, a cavity length sweep is presented at nearly constant values of M e , Re 0 
and H/8. Turbulent bump factor image data are presented in Figure 5 and Figure 6. These demonstrate how the flow 
physics can change with cavity length. In contrast to the laminar data (see Ref. 8), significant vortex interactions 
occur on the sides at the top of the turbulent cavity as the flow strongly expands into the cavity, carrying high 
temperature fluid. The images show the significant increase in floor heating on the center and in the corner where 
vortex scrubbing elevates the temperatures. The complexity of the flow field is further evidenced by the occurrence 
of multiple heating peaks on the downstream endwall under some conditions, for example see Figure 5. In fact, the 
peak heating may not occur on the centerline because of this flow characteristic. 

Linecut plots of the BF distributions are presented in each of the image figures for illustration purposes only. 
These are taken longitudinally along the cavity centerline, spanwise at the cavity entrance, spanwise near the mid- 
cavity location, and spanwise at the cavity exit. The spanwise spread and the longitudinal decay of the cavity wake 
are also highlighted as flow recovers to an undisturbed state. 

V. Summary of Laminar Bump Factor Correlations 

Previously developed laminar correlations are presented here in preparation for the turbulent development. 
Everhart 8 graphically demonstrated that neither BF AVG nor BF M ax could be correlated with the individual primary 
test variables (e.g. M e , Re 0 , Re 0 /M e , L/H, H/8, L/8) for laminar flow. Accordingly, the test variables were combined 
into laminar correlation parameters to describe the variation of heating over a wide range of conditions for both open 
and closed cavity flow physics using aggregate and separate Mach 6 and Mach 10 data. These correlation 
parameters are given by 
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X = In 


M° e Re T e 





( 1 ) 


For later comparison with the turbulent results, the laminar Mach 6 power coefficients are presented in Table 2. 

The correlation parameters for BF A vg and BF M ax were developed using the Mach 6 baseline, rectangular cavity 
data from T6888. These parameters were then input to the commercially available SYSTAT TableCurve 2D program 
to establish the mathematical relationships for the heating variations. The coefficients a and b in the linear (in 
transformed log-log space) correlation curve 


Y = a + bX (2) 

were determined by the program. Output from the program included Confidence Limits on the curve coefficients 
and the Linear Correlation Coefficient, R 2 , as statistical measures of the goodness-of-fit to the data. Confidence 
Limits are a measure of the uncertainty in the average value of a coefficient drawn from an existing population. The 
Linear Correlation Coefficient measures how well the data fit a linear mathematical model, with R 2 =l providing a 
perfect fit. Prediction Limits measure the ability of the curve to estimate the value of a future observation obtained 
from the population. By statistical necessity, Prediction Limits must be wider than Confidence Limits to capture the 
additional uncertainty carried by the new data. A discussion of these statistical concepts is provided in Morrison 28 . 

The Mach 6 laminar correlation curves and the Mach 6 data are presented in Figure 7 for reference. Open cavity 
data are represented by open symbols and the closed cavity data are represented by closed symbols. All of the data 
scatter uniformly about the correlation curve within the bounds of the 99% Prediction Limits. Because these data 
were acquired on a laminar low- temperature flat plate, the scatter band is wider than would be anticipated for a 
higher temperature turbulent flow. Therefore, it is anticipated that the prediction intervals for the turbulent 
correlations will be narrower. For clarity, only the laminar curves will be presented for comparison in subsequent 
figures. 


VI. Turbulent Closed-Cavity Heating Correlations 

Beginning with the laminar relationships presented in the previous section, the first step in analyzing the 
turbulent data was to assess whether Eq. 1 and parameter powers in Table 2 yield any significant correlation of 
BFavg and BF M ax- Direct, unadjusted application of these expressions is provided in Figure 8, revealing well- 
defined offset variations, but with significant scatter that is functionally related to the cavity depth. Therefore, the 
values of the laminar powers provided in Table 2 must be adjusted to remove the depth effect for the turbulent 
conditions. 

New turbulent power coefficients were developed following the process developed in Ref. 8. The laminar 
powers were iteratively changed to maximize the value of the Correlation Coefficient, R 2 . These new turbulent 
values are presented in Table 4 and the results of the modifications are shown in Figure 9. Significant improvements 
in the linearity are provided with the adjusted values that are indicative of the reduced scatter provided by the 
correlation. The R 2 for BF AVG has increased from 0.950 to 0.993; the R 2 for BF M ax has increased from 0.930 to 
0.988. The Degrees of Freedom Adjusted R 2 values are 0.993 for BF AVG and 0.987 for BF M ax, all measures 
implying a high degree of linear modeling. 

Convergence to a local minimum/maximum is a concern in all optimization problems. A correlation 
development question was whether a single set of correlation parameters could be developed to describe both 
laminar and turbulent conditions. As a check on the minimization process, the laminar data were tested against the 
turbulent correlations. The turbulent parameters did describe the gross behavior of the laminar heating, but did so 
poorly. Therefore, separate closed-cavity correlations of the heating are required for the two boundary layer states. 

As before, correlation curve coefficient values were developed and they are provided in Table 5. Also provided 
are the 99% Upper Confidence Limits on the curve fit, and the 99% Upper Prediction Limits. Interestingly, the 
turbulent correlations are nearly parallel to the laminar correlations, as determined by both visual observation in 
Figure 9 and by comparison of the slopes (b coefficient) given in Table 3 with those in Table 5. The turbulent 99% 
uncertainty bounds are narrower than those developed for the laminar correlations because the higher heating 
reduces the uncertainty in the individual measurements. Fewer data are available over a narrower test space resulting 
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in a broadening of the uncertainty at the extents of the turbulent correlation curve. Therefore, the uncertainty was 
estimated using a linear approximation between the end points. 

Final expressions for the correlations are obtained by substituting equations (1) into Eq. (2), yielding 



BF = e a 


M° e R e z e 



(3) 


where the values of the powers a, (3, a, and t are obtained from Table 4, and the a and b coefficients are obtained 
from Table 5. 

A relationship for the Bump Factor uncertainty, ABF/BF, can be obtained by subtracting the linear relationship 
for the Correlation Curve from the linear relationship from the Upper Prediction Limit. This is given by 


ABF 

BF 


— g( a UPL ~ a ) 


M G e Re; 



-1 


(4) 


where ABF=BF UPL -BF and, as above, the coefficients are obtained from Table 5. Uncertainties in the turbulent 
closed-cavity correlations are given as 


( ABF 

) =1.3457 

Mf 5 Re°- 35 fU 

l BF , 

' AVGJurb 

L UyJ 


-1 = 0.35 


and 


( 5 ) 


fABF 

= 1.6991 

M"Rer^l 

l BF y 

‘ MAX, Turb 

L UJJ 


( 6 ) 


Reflecting the lower temperature measurements, larger uncertainties were present for the laminar correlations 
developed in Ref. 8. For comparison, they are 


ABF 

BF 


AVG,Lam 


= 0.75 


(V 


and 


= 1.28 

MAX, Lam ( 8 ) 

VII. Turbulent Open-Cavity Heating Correlations 

Fully turbulent closed-cavity heating data were not discovered in the open literature; however, several sources of 
turbulent open-cavity data with sufficient analysis information for obtaining cavity- floor averages were available. 
Because all of the archival data were acquired using discrete sensors, measurements in the vicinity of the endwall 
top edge were difficult to obtain, allowing only extrapolated peak values containing large uncertainties. 
Accordingly, only correlations using the floor averages (an integral process with less inherent error) were attempted. 
The data sources where the boundary layer entering the cavity is turbulent are Charwat 29,30 , Emery 31 , Netterfield 32 , 
and Thomann 33 . These data were acquired on the wind tunnel floor and on an axisymmetric cone. Data acquired by 
Hahn 34 (axisymmetric cylinder) and Nestler 35 (flat plate) have transitional boundary layers entering the cavity. 
According to Nestler, the heating augmentation for transitional flow will be the same as that for turbulent flow over 
a separated cavity. This important hypothesis for cavity design and impact damage assessment will be demonstrated 
as erroneous. 

These archival data represent an extremely wide range of turbulent flow and geometric conditions that are 
tabulated in Table 6 and plotted as BF A vg versus L/8 in Figure 10. Cavity-entry Mach number varies from 1.8 to 9; 
momentum thickness Reynolds number varies from 2098 to 32288; depth H/8 varies from 0.44 to 4.17; and length 


ABF 

BF 
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L/8 varies from 0.33 to 44. All of the cavities have a rectangular (or nearly rectangular) profile geometry with a 
width W/H that varies from 1.25 for the three dimensional cases to infinity for the axisymmetric cases. In many 
cases the required boundary layer parameters were unavailable in the reports. Referencing White 36 , estimates of the 
turbulent thicknesses were obtained from 


- = 0.37Re: 1/5 and - = 0.036Re: 1,s 

x s x s (9) 

where Re s is the Reynolds number based on boundary layer run length. Combining Eqs. (9) yields 


| = 0.0 9 7 


( 10 ) 


Equation 10 was used to determine the momentum thickness Reynolds number where required. Average floor- 
heating data was obtained by integrating cavity-heating distributions when not provided directly in report figures. 

An initial attempt at correlating the turbulent open cavity BF A vg using the turbulent closed cavity correlation is 
shown in Figure 11. These comparisons were developed using Eqs. (1), and powers from Table 4. While the 
turbulent closed-cavity correlation collapses each individual data set, the different sets are unacceptably separated. 
Unlike the single-curve correlation provided for laminar open and closed cavities in Ref. 8, two distinct correlations 
are required for turbulent open and closed cavities. A successful direct application of the laminar correlation powers 
(Table 2) to the turbulent open cavity data is presented in Figure 12. No attempt at further scatter-reduction via 
modification of the laminar powers was attempted. The implication of this collapse is that cavity heating is a vortex 
and gradient driven problem where the flow gradients for laminar open and closed cavities are relatively benign, as 
they are for open turbulent cavities. Thus, the flows scale similarly, providing a single set of correlation powers. A 
definitive statement regarding the behavior of BF M ax cannot be made without appropriate data; however, similar 
behavior is anticipated. Closed turbulent cavities provide three distinct internal flow fields that develop severe 
turning gradients. These differences result in much higher heating levels and greater variability along the cavity 
floor, and they are reflected in the powers in the scaling parameters and in the turbulent curve-fit slopes. With the 
correct turbulent scaling defined, it is now observationally clear that transitional-entry cavities are outliers, and the 
data should be excluded. 

Finally, Figure 13 presents the turbulent open-cavity correlation curves. The powers and curve coefficients are 
presented in Table 7 and Table 8, respectively. As with the turbulent closed cavity, curve-fit coefficients, 
Confidence Fimits, and Prediction Fimits were determined at the 99% level of uncertainty using the SYSTAT 
TableCurve 2D program. Analytical expressions for the open-cavity uncertainties can be obtained following the 
procedure presented in Eqns. (3) through (6). 


VI 1 1 . Concluding Remarks 

A massive cavity-heating database was generated to support tool development for impact damage assessment of 
re-entry vehicle thermal protection systems. These data were previously documented and they formed the basis of 
the current laminar analysis methods that are incorporated into the Cavity Heating Tool (Version 3). These laminar 
cavity data were recently re-analyzed and new correlations with quantifiable uncertainty levels were developed and 
reported for the floor-averaged bump factor (or heating augmentation) and the maximum endwall bump factor. 

Included in the dataset were turbulent cases for comparison with laminar cavity heating and for evaluation of the 
existing turbulent analysis methodologies. These turbulent data have now been reduced using new turbulent 
simulations of the base model (no cavity) flow field. Application of the laminar-cavity heating correlations to the 
turbulent closed-cavity data captured the variations of the cavity-floor-averaged heating and the maximum heating 
on the cavity endwall; however, data scatter was larger than desired. New turbulent coefficients were successfully 
derived, and they are very descriptive of cavity heating under common surface flow and geometric variations. 
Archival turbulent open-cavity data covering an extremely wide range of flow conditions and geometry have been 
retrieved from the literature and they were successfully correlated. Current and future vehicle programs should 
benefit from their application to impact-damage assessment and to vehicle risk reduction efforts during design. 
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Tables 


Table 1.- Turbulent test matrix for 20-Inch Mach 6 Tunnel Test 6888. 


s 

3 

Pi 

Config 

a ( deg ) 

Re 

L 

( ft 1 ) 

* 'S 

" fl 
o .2 

x ^ 

M e 

Re 0 

Ee e 

M e 

8 

( inch ) 

H 

8 

w 

H 

L 

H 

L 

8 

0 

> 

< 

to 

CP 

< 

to 

CP 

39 

27 b 

0 

2.0 

10.5 

2.93 

123 

42 

0.248 

— 

— 

— 

— 

— 

— 

31 

27 b 

0 

4.0 

18 

2.93 

150 

50 

0.312 

— 

— 

— 

— 

— 

— 

42 

27 b 

-10 

2.0 

10.5 

2.31 

155 

67 

0.178 

— 

— 

— 

— 

— 

— 

226 

— 

-5 

0.6 

10.5 

2.63 

97 

37 

0.299 

— 

— 

— 

— 

— 

— 

215 

— 

-5 

1.9 

10.5 

2.62 

135 

51 

0.216 

— 

— 

— 

— 

— 

— 

117 

28 

0 

2.0 

10.5 

2.93 

123 

42 

0.248 

0.585 

2.6 

15.6 

9.1 

1.143 

5.454 

95 

29 

0 

2.0 

10.5 

2.93 

123 

42 

0.248 

0.612 

2.5 

18.1 

11.1 

0.895 

5.084 

67 

30 

0 

2.0 

10.5 

2.93 

123 

42 

0.248 

0.626 

2.5 

20.8 

13.0 

1.007 

5.479 

87 

31 

0 

2.0 

10.5 

2.93 

123 

42 

0.248 

0.618 

2.5 

26.5 

16.4 

0.945 

4.668 

108 

32 

0 

2.0 

10.5 

2.93 

123 

42 

0.248 

0.578 

2.6 

34.0 

19.7 

0.791 

4.939 

ID 

36 

0 

2.0 

10.5 

2.93 

123 

42 

0.248 

1.389 

2.4 

14.3 

19.8 

1.264 

8.425 

m 

37 

0 

2.0 

10.5 

2.93 

123 

42 

0.248 

1.321 

2.6 

18.2 

24.1 

1.249 

5.741 

mm 

38 

0 

2.0 

10.5 

2.93 

123 

42 

0.248 

1.386 

2.4 

20.5 

28.4 

1.188 

8.270 


39 

0 

2.0 

10.5 

2.93 

123 

42 

0.248 

1.354 

2.5 

26.2 

35.5 

0.999 

8.446 

■a 

40 

0 

2.0 

10.5 

2.93 

123 

42 

0.248 

1.404 

2.4 

30.5 

42.8 

0.923 

6.378 

ID 

41 

0 

4.0 

18 

2.99 

150 

50 

0.312 

0.535 

2.1 

12.8 

6.8 

1.311 

5.585 


42 

0 

4.0 

18 

2.99 

150 

50 

0.312 

0.458 

2.5 

17.9 

8.2 

1.325 

4.754 

mm 

43 

0 

4.0 

18 

2.99 

150 

50 

0.312 

0.474 

2.4 

20.6 

9.8 

1.380 

4.550 

D9 

44 

0 

4.0 

18 

2.99 

150 

50 

0.312 

0.459 

2.5 

26.5 

12.1 

1.267 

5.990 

ID 

48 

0 

4.0 

18 

2.99 

150 

50 

0.312 

1.069 

2.4 

13.8 

14.8 

1.710 

7.509 


49 

0 

4.0 

18 

2.99 

150 

50 

0.312 

1.106 

2.3 

16.3 

18.0 

1.410 

6.019 

ID 

50 

0 

4.0 

18 

2.99 

150 

50 

0.312 

1.111 

2.3 

19.1 

21.2 

1.421 

10.956 

214 

28 

-5 

1.9 

10.5 

2.62 

135 

51 

0.216 

0.671 

2.6 

15.6 

10.5 

1.110 

4.595 

215 

29 

-5 

1.9 

10.5 

2.62 

135 

51 

0.216 

0.701 

2.5 

18.1 

12.7 

1.058 

3.727 

219 

31 

-5 

1.9 

10.5 

2.62 

135 

51 

0.216 

0.709 

2.5 

26.5 

18.8 

1.047 

3.370 

222 

32 

-5 

1.9 

10.5 

2.62 

135 

51 

0.216 

0.663 

2.6 

34.0 

22.5 

1.000 

3.726 

116 

28 

-10 

2.0 

10.5 

2.31 

155 

67 

0.178 

0.814 

2.6 

15.6 

12.7 

1.210 

4.130 

76 

72 

-10 

2.0 

10.5 

2.31 

155 

67 

0.178 

0.813 

7.9 

20.0 

16.3 

0.808 

2.449 

94 

29 

-10 

2.0 

10.5 

2.31 

155 

67 

0.178 

0.851 

2.5 

18.1 

15.4 

1.057 

3.778 

66 

30 

-10 

2.0 

10.5 

2.31 

155 

67 

0.178 

0.871 

2.5 

20.8 

18.1 

1.127 

4.780 

85 

31 

-10 

2.0 

10.5 

2.31 

155 

67 

0.178 

0.860 

2.5 

26.5 

22.8 

1.038 

3.608 

107 

32 

-10 

2.0 

10.5 

2.31 

155 

67 

0.178 

0.805 

2.6 

34.0 

27.3 

1.008 

4.248 

226 

30 

-5 

0.6 

10.5 

2.65 

99 

37 

0.299 

0.518 

2.4 

20.8 

10.8 

0.917 

2.727 
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Table 2 - Laminar cavity correlation coefficient powers derived from Mach 6 baseline data. 


Power 

BFavg 

BF M ax 

a 

3.30 

3.30 

3 

2.70 

3.70 

a 

0.00 

0.20 

T 

0.10 

0.05 


Table 3.- Laminar cavity correlation coefficient values for y=a+bx fit of Mach 6 baseline data. 



BFavg 

BFmax 

a 

b 

DF Adj R 2 

a 

b 

DF Adj R 2 

Correlation 

-2.8251 

2.9650 

0.9882 

-0.6305 

3.8695 

0.9895 

99% UCL 

-2.1834 

2.9715 


0.1707 

3.8737 


99% UPL 

-2.7237 

2.9540 


-0.5137 

3.8588 



Table 4.- Turbulent cavity correlation coefficient powers. 


Power 

BFavg 

BFmax 

a 

2.65 

3.40 

3 

3.20 

4.20 

a 

0.25 

0.80 

T 

0.35 

0.40 


Table 5.- Turbulent cavity correlation coefficient values for y=a+bx fit of Mach 6 baseline data. 



BFavg 

BFmax 

a 

b 

DF Adj R 2 

a 

b 

DF Adj R 2 

Correlation 

-3.8686 

2.7683 

0.9925 

-8.4946 

3.9942 

0.9879 

99% UCL 

-3.4116 

2.8636 


-7.1230 

4.2394 


99% UPL 

-3.5717 

2.7711 


-7.9645 

4.0036 
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Table 6.- Archival cavity-floor heating data, geometry, and conditions. 


Source 

Configuration 

Me 

Red 

8 

(inch) 

H 

(inch) 

W 

(inch) 

L 

(inch) 

H/8 

W/H 

L/H 

L/5 

BFavg 

Charwat 

Tunnel Floor 

2.1 

7943 

0.425 

0.250 

3.00 

0.750 

0.59 

12.00 

3.0 

1.76 

1.100 

Charwat 

Tunnel Floor 

2.9 

6310 

0.540 

0.250 

3.00 

0.750 

0.46 

12.00 

3.0 

1.39 

0.900 

Charwat 

Tunnel Floor 

3.5 

4409 

0.568 

0.250 

3.00 

0.750 

0.44 

12.00 

3.0 

1.32 

0.850 

Charwat 

Tunnel Floor 

2.1 

7943 

0.425 

0.250 

3.00 

1.500 

0.59 

12.00 

6.0 

3.53 

1.090 

Charwat 

Tunnel Floor 

2.9 

6310 

0.540 

0.250 

3.00 

1.500 

0.46 

12.00 

6.0 

2.78 

0.920 

Charwat 

Tunnel Floor 

3.5 

4409 

0.568 

0.250 

3.00 

1.500 

0.44 

12.00 

6.0 

2.64 

0.860 

Charwat 

Tunnel Floor 

2.1 

7943 

0.425 

0.250 

3.00 

2.250 

0.59 

12.00 

9.0 

5.29 

1.150 

Charwat 

Tunnel Floor 

2.9 

6310 

0.540 

0.250 

3.00 

2.250 

0.46 

12.00 

9.0 

4.17 

0.980 

Charwat 

Tunnel Floor 

3.5 

4409 

0.568 

0.250 

3.00 

2.250 

0.44 

12.00 

9.0 

3.96 

0.880 

Charwat 

Tunnel Floor 

2.1 

7943 

0.425 

0.500 

3.00 

1.500 

1.18 

6.00 

3.0 

3.53 

0.680 

Charwat 

Tunnel Floor 

2.9 

6310 

0.540 

0.500 

3.00 

1.500 

0.93 

6.00 

3.0 

2.78 

0.650 

Charwat 

Tunnel Floor 

3.5 

4409 

0.568 

0.500 

3.00 

1.500 

0.88 

6.00 

3.0 

2.64 

0.640 

Charwat 

Tunnel Floor 

2.1 

7943 

0.425 

0.500 

3.00 

3.000 

1.18 

6.00 

6.0 

7.06 

0.660 

Charwat 

Tunnel Floor 

2.9 

6310 

0.540 

0.500 

3.00 

3.000 

0.93 

6.00 

6.0 

5.56 

0.530 

Charwat 

Tunnel Floor 

3.5 

4409 

0.568 

0.500 

3.00 

3.000 

0.88 

6.00 

6.0 

5.29 

0.500 

Charwat 

Tunnel Floor 

2.1 

7943 

0.425 

0.500 

3.00 

4.500 

1.18 

6.00 

9.0 

10.59 

0.580 

Charwat 

Tunnel Floor 

2.9 

6310 

0.540 

0.500 

3.00 

4.500 

0.93 

6.00 

9.0 

8.33 

0.520 

Charwat 

Tunnel Floor 

3.5 

4409 

0.568 

0.500 

3.00 

4.500 

0.88 

6.00 

9.0 

7.93 

0.420 

Charwat 

Tunnel Floor 

2.1 

7943 

0.425 

0.250 

3.00 

3.000 

0.59 

12.00 

12.0 

7.06 

1.250 

Charwat 

Tunnel Floor 

2.9 

6310 

0.540 

0.250 

3.00 

3.000 

0.46 

12.00 

12.0 

5.56 

1.120 

Charwat 

Tunnel Floor 

3.5 

4409 

0.568 

0.250 

3.00 

3.000 

0.44 

12.00 

12.0 

5.29 

0.940 

Emery 

Tunnel Floor 

2.9 

32288 

0.150 

0.250 

0.75 

0.75 

1.67 

3.00 

3.0 

5.00 

0.602 

Emery 

Tunnel Floor 

2.9 

32288 

0.150 

0.250 

0.75 

1.00 

1.67 

3.00 

4.0 

6.67 

0.563 

Emery 

Tunnel Floor 

2.9 

32288 

0.150 

0.250 

0.75 

1.25 

1.67 

3.00 

5.0 

8.33 

0.462 

Emery 

Tunnel Floor 

2.9 

32288 

0.150 

0.250 

0.75 

1.50 

1.67 

3.00 

6.0 

10.00 

0.434 

Emery 

Tunnel Floor 

2.9 

32288 

0.150 

0.250 

0.75 

1.75 

1.67 

3.00 

7.0 

11.67 

0.403 

Emery 

Tunnel Floor 

2.9 

32288 

0.150 

0.250 

0.75 

2.00 

1.67 

3.00 

8.0 

13.33 

0.385 

Emery 

Tunnel Floor 

2.9 

32288 

0.150 

0.250 

0.75 

2.25 

1.67 

3.00 

9.0 

15.00 

0.407 

Emery 

Tunnel Floor 

2.9 

32288 

0.150 

0.250 

0.75 

2.50 

1.67 

3.00 

10.0 

16.67 

0.475 

Emery 

Tunnel Floor 

2.9 

32288 

0.150 

0.250 

0.75 

2.75 

1.67 

3.00 

11.0 

18.33 

0.478 

Hahn 

Axisym_Cyl 

6.3 

4188 

4.584 

19.100 


101.80 

4.17 


5.3 

22.21 

1.255 

Hahn 

Axisym_Cyl 

6.3 

2797 

4.584 

19.100 


152.40 

4.17 


8.0 

33.25 

1.470 

Hahn 

Axisym_Cyl 

6.3 

2098 

4.584 

19.100 


203.20 

4.17 


10.6 

44.33 

1.800 

Netterfield 

Axisym_Cone 

9.0 

9925 

0.236 

0.098 


0.079 

0.42 


0.8 

0.33 

0.850 

Netterfield 

Axisym_Cone 

9.0 

9925 

0.236 

0.098 


0.157 

0.42 


1.6 

0.67 

0.658 

Netterfield 

Axisym_Cone 

9.0 

9925 

0.236 

0.098 


0.236 

0.42 


2.4 

1.00 

0.520 

Netterfield 

Axisym_Cone 

9.0 

32109 

0.236 

0.098 


0.079 

0.42 


0.8 

0.33 

0.410 

Netterfield 

Axisym_Cone 

9.0 

32109 

0.236 

0.098 


0.157 

0.42 


1.6 

0.67 

0.430 

Netterfield 

Axisym_Cone 

9.0 

32109 

0.236 

0.098 


0.236 

0.42 


2.4 

1.00 

0.500 

Thomann 

Tunnel Floor 

1.8 

12300 

0.424 

0.787 

3.00 

2.362 

1.86 

3.80 

3.0 

5.57 

0.590 

Thomann 

Tunnel Floor 

1.8 

12300 

0.424 

0.394 

3.00 

2.362 

0.93 

7.60 

6.0 

5.57 

0.554 

Thomann 

Tunnel Floor 

1.8 

12300 

0.424 

0.197 

3.00 

2.362 

0.46 

15.20 

12.0 

5.57 

0.709 

Nestler 

Flat Plate 

8.5 

9571 

0.320 

1.000 

4.00 

5.0 

4.17 

3.13 

5.0 

15.6 

0.515 

Nestler 

Flat Plate 

8.5 

9571 

0.320 

0.800 

4.00 

12.0 

4.17 

2.50 

15.0 

37.5 

1.211 

Nestler 

Flat Plate 

8.5 

9571 

0.320 

0.400 

4.00 

12.0 

4.17 

1.25 

30.0 

37.5 

0.702 


12 







Table 7.- Turbulent open cavity correlation coefficient powers. 


Power 

BF A vg 

a 

3.30 

3 

2.70 

a 

0.00 

T 

0.10 


Table 8.- Turbulent open-cavity correlation coefficient values for y=a+bx fit of archival data. 



BFavg 

. DF Adi 

a b c R 2 

Correlation 
99% UCL 
99% UPL 

-2.3741 2.6034 0.0000 0.9752 

-1.8950 2.3604 0.0538 

-1.0967 2.5399 0.0139 


Figures 





Figure 3. Sample cavity insert model used in Mach 6 for 
T6888. 
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Heating Distribution 



(a) Computational and experimental heating 
distributions. 



(b) Surface flow properties. 


Figure 4.- Turbulent CFD simulations for Test 6888 Run 31. 


Test 6888/1 48 Configuration 49 
L/H = 1 6.3 W/H = 2.3 H/8 = 1.106 
M e = 2.99 Re 0 = 1 50 

Endwall 



Figure 5.- Effect of length variation for M e =2.99, Re e =150, H/8=l.l — L/H~16. 
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Test 6888/147 Configuration 50 
L/H = 19.1 W/H = 2.3 H/8 = 1.111 



Figure 6.- Effect of length variation for M e =2.99, Re e =150, and H/8=l.l — 1/11=19. 




a) Cavity floor correlation. b) Cavity endwall maximum correlation. 

Figure 7.- Laminar baseline Mach 6 correlation. 
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Figure 8.- Correlation of turbulent data with laminar correlation coefficients. 
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Figure 9.- Correlation of turbulent data with turbulent correlation coefficients, including 99% uncertainty 

analysis. 
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Figure 10.- Archival turbulent-entry open-cavity floor-averaged heating augmentation data. 



Figure 11.- Turbulent-entry open-cavity floor-heating data scaled using turbulent closed-cavity correlation. 
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Figure 12.- Correlated turbulent cavity-floor-averaged heating augmentation. 
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Figure 13.- Turbulent open-cavity floor-averaged heating augmentation. 
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